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Abstract 

Jet substructure and differential cross sections for jets produced in the photo- 
production and deep inelastic ep scattering regimes have been measured with 
the ZEUS detector at HERA using an integrated luminosity of 82.2 pb^^. 
The substructure of jets has been studied in terms of the jet shape and sub- 
jet multiplicity for jets with transverse energies ^j?* > 17 GeV. The data 
are well described by the QCD calculations. The jet shape and subjet mul- 
tiplicity are used to tag gluon- and quark-initiated jets. Jet cross sections 
as functions of -E;^*, jet pseudorapidity, the jet-jet scattering angle, dijet in- 
variant mass and the fraction of the photon energy carried by the dijet sys- 
tem are presented for gluon- and quark-tagged jets. The data exhibit the be- 
haviour expected from the underlying parton dynamics. A value of as{Mz) of 
as{Mz) = 0.1176 ± 0.0009 (stat.) to!^Xl (exp.) IomyI (th.) was extracted from 
the measurements of jet shapes in deep inelastic scattering. 
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1 Introduction 



Jet production in ep collisions provides a fruitful testing ground of perturbative QCD 
(pQCD). Measurements of differential cross sections for jet production [1-9] have allowed 
detailed studies of parton dynamics, tests of the proton and photon parton distribution 
functions (PDFs) as well as precise determinations of the strong coupling constant, a^. 
Most of these measurements refer to the production of jets irrespective of their partonic 
origin - quarks or gluons - and, therefore, have only provided general tests of the partonic 
structure of the short- distance process and of combinations of the proton and/or photon 
PDFs. The identification of quark- and gluon-initiated jets would allow more stringent 
tests of the QCD predictions. Such measurements of the production of jets containing 
a heavy quark have been made by means of tagging specific decay channels [10]. In 
the present study, quark- and gluon-initiated jets are identified on a statistical basis by 
utilising their internal structure. 

Two kinematic regimes have been studied: photoproduction (7p) and neutral current 
(NC) deep inelastic ep scattering (DIS). Photoproduction at HERA is studied by means 
of ep scattering at low four- momentum transfers {Q"^ ~ 0, where is the virtuality of 
the exchanged photon). In photoproduction, two types of QCD processes contribute to 
jet production at leading order (LO) [11, 12]: either the photon interacts directly with a 
parton in the proton (the direct process) or the photon acts as a source of partons which 
scatter off those in the proton (the resolved process). Jet production in NC DIS up to 
LO in as proceeds as in the quark-parton model {Vq — > q, where V = 'y or Z^) or via the 
boson-gluon fusion (Vg — > qq) and QCD-Compton {Vq — > qg) processes. 

This paper is organised as follows. Section |21 gives the theoretical expectations for the 
measurements presented. The experimental set-up and data selection are described in Sec- 
tions El and El respectively. Section explains the QCD calculations used in this analysis. 
The corrections applied to the data and systematic uncertainties are given in Section IHl 
The results on the mean integrated jet shape and subjet multiplicity in photoproduction 
and NC DIS are presented in Section [71 The measurements of differential inclusive jet 
cross sections as a function of the jet pseudorapidity^ , rj^^^, and jet transverse energy, E^r^^, 
for samples of jets in the photoproduction and NC DIS regimes, separated according to 
their shape and subjet multiplicity, are presented in Section |H1 In addition, measurements 
of differential dijet cross sections in photoproduction as a function of cos 6*, where 6* is the 
angle between the jet-jet axis and the beam direction in the dijet centre-of-mass system, 
the dijet invariant mass, and the fraction of the photon momentum participating in 

^ The ZEUS coordinate system is a right-handed Cartesian system, with the Z axis pointing in the 
proton beam direction, referred to as the "forward direction" , and the X axis pointing left towards 
the centre of HERA. The coordinate origin is at the nominal interaction point. 
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the production of the two jets with highest -E;^*, x^^, are also presented for a variety of 
tagged-jet configurations. The results are compared to leading-logarithm parton-shower 
calculations and used to investigate the dynamics underlying the production of specific 
tagged-jet final states. Finally, in Section IHl the measurements of the mean integrated 
jet shape in NC DIS are compared to next-to-leading-order (NLO) QCD predictions and 
used to extract a^. 

2 Theoretical expectations 

The internal structure of a jet depends mainly on the type of primary parton - quark or 
gluon - from which it originated and to a lesser extent on the particular hard scattering 
process. At sufficiently high jet transverse energy, where the influence of fragmentation 
effects becomes negligible, the internal structure of a jet is calculable in pQCD. Such 
calculations predict that gluon-initiated jets are broader than quark-initiated jets due to 
the larger colour charge of the gluon. The jet shape [13] and subjet multiplicity [14] 
can be used to study the internal structure of the jets and to classify them: a "broad" - 
jet sample is enriched in gluon-initiated jets, whereas a "narrow"-jet sample is enriched 
in quark-initiated jets. Thus, measurements of cross sections for broad- and narrow-jet 
samples allow the contributing hard-scattering subprocesses to be disentangled. 

The dominant partonic subprocesses responsible for jet photoproduction in the kinematic 
region presented in this paper are 7(7 —>■ qq and q^gp qg, where q^ (gp) denotes a 
quark (gluon) coming from the photon (proton). The kinematics of these two-to-two 
subprocesses are such that the majority of the jets in the region rf'^^ < originate from 
outgoing quarks, whereas the fraction of gluon-initiated jets increases as rf'^^ increases. 

The distribution in 6* reflects the underlying parton dynamics and is sensitive to the 
spin of the exchanged particle. In the case of direct-photon interactions, the contributing 
subprocesses at LO QCD are (i) 7g(g) gq{q) and (ii) •yg — ^ qq, which involve quark 
exchange. The behaviour of the dijet angular distribution as | cos6'*| — ^ 1 is the same for 
all direct subprocesses and proportional to (1 — | cos 6*\)^^. In the case of resolved-photon 
interactions, the contributing subprocesses are qg qg, qq' qq', gg gg, ... . The 
dominant subprocesses are those that involve gluon exchange and the behaviour of the 
dijet angular distribution as | cos6'*| — ^ 1 is proportional to (1 — | cos^*|)~^. The different 
behaviour of the dijet angular distribution for resolved and direct processes has been 
measured in photoproduction at HERA [3,15]. The study of the angular distribution for 
dijet events with tagged quark- and/or gluon-initiated jets in the final state, provides then 
a handle to investigate the underlying parton dynamics further. 

Measurements of jet substructure in NC DIS allow a determination of a^. In zeroth-order 



2 



pQCD, a jet consists of only one parton and the jets have no substructure. The first 
non-trivial contribution to the jet substructure is given by 0{as) processes in which, e.g., 
a quark radiates a gluon at a small angle; these are proportional to the rate of parton 
emission and thus to as- For DIS in the laboratory frame, all necessary QCD corrections 
to the jet cross sections for the determination of from the jet substructure are available. 

2.1 Jet-shape and sub jet-multiplicity definitions 

The /ct cluster algorithm [16] was used in the longitudinally invariant inclusive mode [17] 
to define jets in the hadronic final state. The integrated jet shape, ip{r), is defined using 
only those particles belonging to the jet as the fraction of the jet transverse energy that 
lies inside a cone in the rj — (f plane of radius r concentric with the jet axis: 



where ET{r) is the transverse energy within the given cone of radius r. The mean inte- 
grated jet shape, {ip{r)), is defined as the averaged fraction of the jet transverse energy 
inside the cone r: 



where A'jets is the total number of jets in the sample. 

The integrated jet shape is calculated at LO in pQCD as the fraction of the jet transverse 
energy, due to parton emission, that lies in the cone segment between r and R = 1: 



dictions of the integrated jet shape, the numerator in the above formula is calculated to 
and the denominator to 0{aas). 

Studies of QCD using jet production in NC DIS at HERA are usually performed in the 
Breit frame. The analysis of jet shapes presented here was performed in the laboratory 
frame, since calculations of this observable in the Breit frame can, at present, only be 
performed to 0{as), precluding a reliable determination of as- However, calculations of 
the jet shape can be performed up to 0(a;^) in the laboratory frame. Furthermore, the 
analysis was performed in the kinematic region defined by > 125 GeV^ since, at lower 
values of Q^, the sample of events with at least one jet with i?;^* > 17 GeV is dominated by 
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dijet events. The calculation of the integrated jet shape for dijet events can be performed 
only up to 0{as), which would severely restrict the accuracy of the predictions. 

Subjets were resolved within a jet by considering all particles associated with the jet and 
repeating the application of the /c-r cluster algorithm until, for every pair of particles i 
and j the quantity dij = mm{ET,i, Etj)^ ■ {{jji — rjjY + {^pi — ^pj)^), where -Et.i, Vi 
ipi are the transverse energy, pseudorapidity and azimuth of particle i, respectively, was 
greater than dcut = UcntiE'^^^Y ■ AH remaining clusters were called subjets. The subjet 
multiplicity, ?T.subjet; depends upon the value chosen for the resolution parameter ycut- 

The 

mean subjet multiplicity, (risubjot), is defined as the average number of subjets contained 
in a jet at a given value of ycut'- 



where 'T'subjet(2/cut) is the number of subjets in jet i. 

3 Experimental set-up 

The data used in this analysis were collected during the 1998-2000 running period, when 
HERA operated with protons of energy Ep = 920 GeV and electrons or positrons^ of 
energy Eg = 27.5 GeV, and correspond to an integrated luminosity of 82.2 ±1.9 pb~^ 

A detailed description of the ZEUS detector can be found elsewhere [18,19]. A brief outline 
of the components that are most relevant for this analysis is given below. Charged particles 
are tracked in the central tracking detector (CTD) [20], which operates in a magnetic 
field of 1.43 T provided by a narrow superconducting solenoid. The CTD consists of 
72 cylindrical drift-chamber layers, organized in nine superlayers covering the polar-angle 
region 15° < ^ < 164°. The transverse- momentum resolution for full-length tracks can 
be parameterised as a{pT)/pT = 0.0058pT © 0.0065 © 0.0014/pr, with pr in GeV. The 
tracking system was used to measure the interaction vertex with a typical resolution along 
(transverse to) the beam direction of 0.4 (0.1) cm and to cross-check the energy scale of 
the calorimeter. 

The high-resolution uranium-scintillator calorimeter (CAL) [21] covers 99.7% of the total 
solid angle and consists of three parts: the forward (FCAL), the barrel (BCAL) and 
the rear (RCAL) calorimeters. Each part is subdivided transversely into towers and 
longitudinally into one electromagnetic section (EMC) and either one (in RCAL) or two (in 

^ Here and in the following, the term "electron" denotes generically both the electron (e^) and the 
positron (e^). 
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BCAL and FCAL) hadronic sections (HAC). The smallest subdivision of the calorimeter 
is called a cell. Under test-beam conditions, the CAL single-particle relative energy 
resolutions were a{E)/E = 0.18/\/E for electrons and a{E)/E = 0.35/v^ for hadrons, 
with E in GeV. 

The luminosity was measured from the rate of the bremsstrahlung process ep — > e'yp. The 
resulting small-angle energetic photons were measured by the luminosity monitor [22], a 
lead-scintillator calorimeter placed in the HERA tunnel a.t Z = —107 m. 

4 Data selection and jet search 

A three-level trigger system was used to select events online [19,23]. At the first level, 
events were triggered by a coincidence of a regional or transverse energy sum in the CAL 
and at least one track from the interaction point measured in the CTD. At the second 
level, a total transverse energy of at least 8 GeV, excluding the energy in the eight CAL 
towers immediately surrounding the forward beampipe, was required, and cuts on CAL 
energies and timing were used to suppress events caused by interactions between the 
proton beam and residual gas in the beampipe. At the third level, a jet algorithm was 
applied to the CAL cells and jets were reconstructed using the energies and positions of 
these cells. Events with at least one (two) jet(s) with Ej- > 10 (6) GeV and rj < 2.5 were 
accepted for the inclusive jet (dijet) samples. For systematic trigger studies, all events 
with a total transverse energy of at least 25 GeV, excluding the energy in the eight CAL 
towers immediately surrounding the forward beampipe, were accepted. No jet algorithm 
was applied in this case. 

In the offline selection, a reconstructed event vertex consistent with the nominal inter- 
action position was required and cuts based on the tracking information were applied to 
reduce beam-induced interactions and cosmic-ray events. The main steps of the selection 
of photoproduction and DIS events are briefly explained below. 

4.1 Selection of the photoproduction sample 

Events from collisions between quasi-real photons and protons were selected offline using 
similar criteria to those reported in a previous publication [2,5]. Charged current DIS 
events were rejected by requiring the total missing transverse momentum, p™'****, to be 
small compared to the total transverse energy, E*?*, p^'^^ ^/E^ < 2 VCeV. Any NC 
DIS events with an identified scattered-electron candidate in the CAL [24] were removed 
from the sample using the method described previously [25]. The remaining background 
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from NC DIS events was estimated by Monte Carlo (MC) techniques to be below 0.3% 
and was neglected. 

The selected sample consisted of events from ep interactions with ^ 1 GeV^ and 
a median ^ 10~^ GeV^. The 7p centre-of-mass energy is given by W.yp = ^/sy, 
where y is the inelasticity variable and y/s is the ep centre-of-mass energy, s = iE^Ep. 
The inelasticity variable was reconstructed using the method of Jacquet-Blondel [26], 
yjB = {E — pz)/2E(,, where E is the total CAL energy and pz is the Z component of the 
energy measured in the CAL cells. The value of y was systematically underestimated by 
~ 20% with an r.m.s. of ~ 10%. This effect, which was due to energy lost in the inactive 
material in front of the CAL and to particles lost in the rear beampipe, was satisfactorily 
reproduced by the MC simulation of the detector. The MC event samples were therefore 
used to correct for this underestimation. The photoproduction sample was restricted to 
142 < W^p < 293 GeV [2,5]. 

4.2 Selection of the NC DIS sample 

Events from NC DIS interactions were selected offline using similar criteria to those re- 
ported in a previous publication [27]. The scattered-electron candidate was identified 
using the pattern of energy deposits in the CAL [24]. The energy, E'^, and polar angle, 9e, 
of the electron candidate were also determined from the CAL measurements. The double- 
angle method [28], which uses 6e and an angle 7 that corresponds, in the quark-parton 
model, to the direction of the scattered quark, was used to reconstruct Q^, Qi^a- The 
angle 7 was reconstructed using the CAL measurements of the hadronic final state. 

An electron candidate of energy E'^ > 10 GeV was required to ensure a high and well 
understood electron-finding efficiency and to suppress background from photoproduction. 
The inelasticity variable as reconstructed from the electron, ye, was required to be below 
0.95. This condition removed events in which fake electron candidates from photoproduc- 
tion background were found in the FCAL. The requirements 38 < {E — pz) < 65 GeV, to 
remove events with large initial-state radiation and to reduce further the photoproduc- 
tion background, and p^^^^/ \/E^ < 3 VCeV, to remove cosmic rays and beam-related 
background, were applied. The kinematic range was restricted to Qda > 125 GeV^. 

4.3 Jet search 

The krp cluster algorithm was used in the longitudinally invariant inclusive mode to re- 
construct jets in the hadronic final state from the energy deposits in the CAL cells. For 
DIS events, the jet algorithm was applied after excluding those cells associated with the 
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scattered-electron candidate. The jet search was performed in the rj — (p plane of the 
laboratory frame. The jet variables were defined according to the Snowmass conven- 
tion [29]. The jets reconstructed from the CAL cell energies are called calorimetric jets 
and the variables associated with them are denoted by -E'^cai' V^cl V^cai- ^ total of 
199237 (98240) events with at least one jet satisfying E^^^^^ > 13 GeV and -1 < r/^^^* < 2.5 
were selected in the photoproduction (DIS) sample. 

5 QCD calculations 

5.1 Leading-logarithm parton-shower Monte Carlo models 

The programs Pythia 6.1 [30] and Herwig 6.1 [31] were used to generate photoproduc- 
tion events for resolved and direct processes. Events were generated using GRV-HO [32] 
for the photon and CTEQ4M [33] for the proton PDFs. In both generators, the par- 
tonic processes are simulated using LO matrix elements, with the inclusion of initial- 
and final-state parton showers. Fragmentation into hadrons is performed using the Lund 
string model [34] as implemented in Jetset [30,35] in the case of Pythia, and a cluster 
model [36] in the case of Herwig. Samples of Pythia including multiparton interactions 
(MI) [37] with a minimum transverse momentum for the secondary scatter of 1 GeV [38] 
were used to study the effects of a possible "underlying event". 

Neutral current DIS events including radiative effects were simulated using the Her- 
acles 4.6.1 [39] program with the Djangoh 1.1 [40] interface to the hadronisation 
programs. HERACLES includes corrections for initial- and final-state radiation, vertex 
and propagator terms, and two-boson exchange. The QCD cascade is simulated using 
the colour-dipole model (CDM) [41] including the LO QCD diagrams as implemented in 
Ariadne 4.08 [42] and, as a systematic check of the final results, with the MEPS model 
of Lepto 6.5 [43]. Both MC programs use the Lund string model for the hadronisation. 
The CTEQ5D [44] proton PDFs were used for these simulations. 

These MC samples were used to correct the data to the hadron level, defined as those 
hadrons with lifetime r > 10 ps. For this purpose, the generated events were passed 
through the ZEUS detector- and trigger-simulation programs based on Geant 3.13 [45]. 
They were reconstructed and analysed by the same program chain as the data. The jet 
search was performed using the energy measured in the CAL cells in the same way as for 
the data. The same jet algorithm was also applied to the final-state particles and to the 
partons available after the parton shower; the jets found in this way are referred to as 
hadronic and partonic jets, respectively. 

Electro weak-radiative, hadronisation and Z'^-exchange effects are not at present included 
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in the NLO QCD programs described in Section l^^ Therefore, samples of MC events were 
generated without electroweak-radiative effects so that the data could be corrected for 
these effects when comparing with the NLO QCD predictions. Additional samples of MC 
events without Z'^-exchange effects were generated to correct the NLO QCD calculations 
for these effects and for hadronisation. 



5.2 NLO QCD calculations 

The NLO QCD calculations of the mean integrated jet shapes in DIS'^ are based on 
the program DiSENT [46]. The calculations use a generalised version of the subtraction 
method [47] and are performed in the massless MS renormalisation and factorisation 
schemes. The number of flavours was set to five; the renormalisation {fin) and factorisation 
(/ii?) scales were both set to fiR = fJ^F = Q', «s was calculated at two loops using = 
220 MeV, which corresponds to as{Mz) = 0.1175. The MRST99 [48] parameterisations of 
the proton PDFs were used as defaults for the comparisons with the data. The calculations 
obtained with DiSENT were cross-checked by using the program Disaster++ [49]. The 
differences were smaller than 0.5% for r > 0.3. 

Since the measurements refer to jets of hadrons, whereas the QCD calculations refer to 
partons, the predictions were corrected to the hadron level using the MC samples described 
in Section Em The multiplicative correction factor, Chad, defined as (^/'(r))had/(^/'('"))par, 
where ('?/'(r))par (had) is the mean integrated jet shape before (after) the hadronisation pro- 
cess, was estimated with both the CDM and MEPS models. The procedure for applying 
hadronisation corrections to the NLO QCD calculations was verified by checking that the 
MC predictions for the integrated jet shape at the parton level reproduced the NLO QCD 
calculations. The agreement was well within 0.2% after adjusting the contributions of the 
different subprocesses (eg —>■ eg, eg — > egg and eg — > eqg) in the MC to reproduce the t]-*^* 
cross section and jet shape of the NLO calculations. The values of Chad obtained with 
the CDM model were taken as the defaults; the predictions from the two models were in 
good agreement. The value of Chad was 0.95 at r = 0.5 for E^r^^ = 21 GeV and approached 
unity as E'^* increased. 



^ Only QCD calculations of jet shapes in DIS are compared to the data because there is no NLO program 
available for similar calculations in photoproduction. 
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6 Corrections and systematic uncertainties 



6.1 Jet energy corrections 

The comparison of the reconstructed jet variables for the hadronic and the calorimetric jets 
in simulated events showed that no correction was needed for rj^^^ and (p^'^^ (r/j"* ~ 77^^^ and 
^p}ct (f'^^^i). However, the transverse energy of the calorimetric jet was an underestimate 
of the corresponding hadronic jet energy by an average of ~ 15%, with an r.m.s. of ~ 10%. 
This underestimation was mainly due to the energy lost by the particles in the inactive 
material in front of the CAL. The transverse-energy corrections to calorimetric jets, as a 
function of ?7^^*[ and -E^^j^j and averaged over (f^^^^, were determined for the photoproduction 
and DIS samples using the corresponding MC-generated events [2,3,27]. Henceforth, jet 
variables without subscript refer to the corrected values. After these corrections to the jet 
transverse energy, events with at least one jet satisfying E^r^^ > 17 GeV and —1 < rj^^^ < 2.5 
were retained for the studies of inclusive jet observables and events with at least two jets 
with E^^" > 17 GeV, E^?'^ > 14 GeV and -1 < r/J^* < 2.5, where the jets are labelled in 
decreasing E^r^^ order, were retained for the dijet studies. 

6.2 Acceptance corrections 

Using the selected data sample of inclusive jets with i?;^* > 17 GeV and —1 < 77-''^* < 
2.5, in the kinematic region defined by < 1 GeV^ and 142 < W^p < 293 GeV in 
photoproduction and > 125 GeV^ in DIS, the mean integrated jet shape and mean 
subjet multiplicity were reconstructed using the CAL cells and corrected to the hadron 
level by MC techniques. The corrected values were determined bin-by-bin using the MC 
samples separately for each region of 77-'®* and i?;^* studied. For this approach to be 
valid, the distributions of the uncorrected integrated jet shape and subjet multiplicity 
in the data must be well described by the MC simulations at the detector level. This 
condition was in general satisfied by the MC photoproduction models. In DIS, to obtain 
the best description of the uncorrected jet shape in the data by the MC simulations, 
the contributions of the different subprocesses (eg — > eg, eg — > egg and eg -h> eqg) were 
reweighted so as to reproduce the mean integrated jet shape and rf'^^ distributions in the 
data. This procedure was applied to the simulations of CDM and MEPS and for each 
region in r]^°^ and E^r^^. The correction factors were evaluated using these tuned versions 
of CDM and MEPS in DIS and the default mixture of resolved and direct processes in 
Pythia, Herwig and Pythia MI. 

Differential jet cross sections were measured using the selected data sample of inclusive 
jet events. Dijet differential cross sections in photoproduction were measured in the same 
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kinematic region as above and refer to the two highest-E';^* jets of hadrons in the event 
with ^j?" > 17 GeV, Ei^^^ > 14 GeV and -1 < r/j^* < 2.5. The Pythia (CDM) MC 
samples of photoproduction (DIS) events were used to compute the acceptance correc- 
tions to the jet distributions. These correction factors took into account the efficiency 
of the trigger, the selection criteria and the purity and efficiency of the jet reconstruc- 
tion. The inclusive jet cross sections were obtained by applying bin-by-bin corrections to 
the measured distributions. The samples of Herwig and MEPS were used to compute 
the systematic uncertainties coming from the fragmentation and parton-shower models in 
photoproduction and DIS, respectively (see Section IU7^ . 



6.3 Systematic uncertainties 

A study of the main sources contributing to the systematic uncertainties of the mea- 
surements was performed. The sources considered in the photoproduction measurements 
are: 

• the effect of the treatment of the parton shower and hadronisation was estimated by 
using the Herwig generator to evaluate the correction factors; 

• the effect of the simulation of the trigger was evaluated by using an alternative trigger 
configuration, as explained in Section EJ in both data and MC events; 

• the effect of the uncertainty on W^p was estimated by varying j/jb by its uncertainty 
of ±1% in simulated events; 

• the effect of the uncertainty on the parameterisations of the proton and photon PDFs 
was estimated by using alternative sets of PDFs in the MC simulation to calculate the 
correction factors; 

• the effect of the uncertainty on the absolute energy scale of the calorimetric jets was 
estimated by varying i?;^* by its uncertainty of ±1% in simulated events. The method 
used was the same as in earlier publications [2,3,50] and verified with the 98-00 data 
sample [5]. 

In the DIS regime, the main sources contributing to the systematic uncertainties of the 
measurements are: 

• the effect of the treatment of the parton shower was estimated by using the MEPS 
model to evaluate the correction factors; 

• the effect of the simulation of the trigger was evaluated by using an alternative trigger 
configuration, as explained in Section |3J in both data and MC events; 
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• the effect of the uncertainty on the scattered-electron identification was estimated by 
using an alternative technique [51] to select the candidates, in both data and MC 
events; 

• the effect of the uncertainty of ±1% in the absolute energy scale of the jets was applied 
to the simulated events; 

• the effect of the uncertainty of ±1% in the absolute energy scale of the scattered- 
electron candidate was applied to the simulated events. 

These uncertainties, for each regime, were added in quadrature to the statistical uncer- 
tainty of the data and are shown as error bars in the figures showing the substructure 
measurements. For the cross-section measurements, the uncertainty arising from that on 
the absolute energy scale of the jets is shown separately. The uncertainty in the luminosity 
determination of 2.25% was not included. 

7 Measurements of jet substructure 
7.1 Jet-shape measurements 

The measured mean integrated jet shape as a function of r, (iplr)), for different regions in 
rf'^^ is shown in Fig.nand Table^for the photoproduction regime. The jets broaden as 77-'*^* 
increases. Leading-logarithm parton-shower predictions from Pythia for resolved plus 
direct processes and gluon- and quark-initiated jets are compared to the data in Fig. [T] 
The measured {ip{r)) is reasonably well described by the MC calculations of Pythia 
for resolved and direct processes for —1 < rf'^^ < 1.5, whereas for 1.5 < rf'^^ < 2.5, 
the measured jets are slightly broader than the predictions. From the comparison with 
the predictions for gluon- and quark-initiated jets, it is seen that the measured jets are 
quark-like for — 1 < tf^^ < and become increasingly more gluon-like as rf'^^ increases. 

Figure 121 shows the same measurements as Fig. Q compared to the predictions of Pythia 
including MI. This model gives rise to jets that are much broader than those observed. The 
predictions using Herwig, also shown in Fig. |2l describe the data well for —1 < rf^^ < 1, 
are slightly narrower than the data for 1 < rf^^ < 1.5 and fail to describe the data for 
1.5 < rf'^^ < 2.5. These results and those presented below are consistent with the previous 
ZEUS study of jet shapes in photoproduction [38] which was performed using an iterative 
cone algorithm and at lower E^r^^. 

Figure El and Tables |21 and 01 show the {ip{r)) in different regions of E^^^ for the photopro- 
duction regime. The jets become narrower as E^r^^ increases. The predictions of Pythia 
for resolved plus direct processes reproduce the data reasonably well. For 17 < -E^* < 29 
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GeV, the predictions for resolved processes alone also describe the data, consistent with 
the dominance of resolved processes in this E}^' region. 

Figure shows the measured mean integrated jet shape at a fixed value of r = 0.5, 
{ijj{r = 0.5)), as a function of rf'^^ in photoproduction. The predictions of Pythia for 
quark-initiated jets lie above the data, while those for gluon-initiated jets lie below the 
data. The prediction of Pythia, including resolved and direct processes, also shown in 
Fig. H^, fails to describe the relatively strong broadening of the measured jet shape for 
7]^'^^ > 1.5. This might be because the fraction of gluon-initiated jets in the region t]^'^^ > 1.5 
is underestimated or that the effects of a possible underlying event in the data have not 
been properly taken into account [38]; however, the prediction of Pythia MI fails to 
describe the data over the whole rf^^ range. Since {ip{r = 0.5)) changes from a value close 
to the upper curve (quark-initiated jets) to a value near the lower curve (gluon-initiated 
jets) as rf^^ increases, the broadening of {tp{r = 0.5)) as rj^^^ increases is consistent with an 
increase of the fraction of gluon-initiated jets. The {ip{r = 0.5)) shows an increase with 
E"^* (see Fig. |3)d). The predictions of Pythia for the dependence of {ip{r = 0.5)) on i?;^* 
in resolved plus direct processes reproduce the data well for i?^* > 21 GeV. Therefore, 
the discrepancies between data and MC are concentrated at low Ei^^ and high rj^'^^ values. 

The measured {ip{r)) for different regions of rj^°^ and E^r^^ is shown in Figs. |S] and Inland 
Tables m to ini for DIS events. Figure [3 shows the measured {ip{r = 0.5)) as a function of 
7]^^^ and E^r^^. There is no significant variation of (^/'(r = 0.5)) with 77^°* in DIS, whereas 
{ijj{r = 0.5)) increases as -E;^* increases, as observed in a previous study [52] using an 
iterative cone algorithm. These conclusions are also in agreement with those of a previous 
publication [27], in which the internal structure of jets in NC DIS was studied using the 
mean subjet multiplicity. 

The NLO QCD calculations of {■ip{r)), corrected for hadronisation and Z°-exchange effects, 
are compared to the data in Figs. to The NLO QCD calculations give a good 
description of (ipir)) for r > 0.2; the fractional differences between the measurements and 
the predictions, also shown in Figs. El and IHl are less than 0.2% for r = 0.5. The sensitivity 
of the measurements to the value of as{Mz) is illustrated in Fig. ^jp by comparing the 
measured {ip^r = 0.5)) to NLO QCD calculations using three different values of as{Mz)- 
The NLO QCD calculations provide a good description of the measured {ip{r = 0.5)) and 
thus this observable was used to determine as{Mz), as explained in Sectional 

Figure IHl shows the rj^°^ and Ei^^ dependence of ('?/'(r = 0.5)) for photoproduction and DIS 
events; the MC predictions of CDM and Pythia for quark- and gluon-initiated jets are 
compared to the data. Figure |Hk shows that the sample of jets in DIS is consistent with 
being dominated by quark-initiated jets with an approximately constant fraction over the 
^jet j^ggiQ^ measured. Photoproduced jets in the backward region are similar to jets in DIS 
and this agreement confirms, independently of the comparison to MC predictions, that 
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they are dominated by quark-initiated jets. The increasing deviation in the integrated jet 
shape for photoproduced jets from that of jets in DIS as 77-'°* (-E't*) increases (decreases) 
can be attributed to the increasing fraction of gluon-initiated jets arising from resolved 
processes. 

7.2 Sub jet-multiplicity measurements 

The measured mean subjet multiphcity as a function of t/cut for different regions of rj^'^^ and 
E^r^^ for photoproduction is shown in Figs. IHl and ^1 and Tables [7|to|ni Figure ITT] shows the 
measured mean subjet multiplicity at a fixed value of ?/cut = 10"^, (nsubjct(l/cut = 10"^)), 
as a function of rj^°^ and E^^^. The measured mean subjet multiplicity increases as rj^^^ 
increases and decreases as E^r^^ increases. 

The predictions of Pythia for quark-initiated, gluon-initiated and all jets at the hadron 
level are compared to the measurements in Figs. IHl to [TTJ The predicted (nsubjet(l/cut)) is 
larger for gluon-initiated jets than for quark-initiated jets in each region of rf^^. For quark- 
or gluon-initiated jets alone, (risubjet(z/cut = 10~^)) exhibits only a small dependence on 
rj^^^ (see Fig. \TTk). The ?7-'*^*-dependence of (nsubjct(ycut = 10~^)) in the calculation for all 
jets is dictated by the rj^^^ variation of the fractions of quark- and gluon-initiated jets. 
This variation, in turn, originates from the different dominant two-body subprocess. The 
calculations using Pythia based on the predicted admixture of quark- and gluon-initiated 
jets give a good description of the measured (r2subjct(2/cut)) as a function of ?/cut5 V'^* and 
E';^*. These results are in agreement with those from the mean integrated jet shape (see 
Section rrTj) . 

8 Study of quark and gluon dynamics 

The predictions of the MC for the jet shape and subjet multiplicity generally reproduce 
the data well and show the expected differences for quark- and gluon-initiated jets. These 
differences are used now to select samples enriched in quark- and gluon-initiated jets to 
study the dynamics of the hard subprocesses in more detail. 

8.1 Selection of quark- and gluon-initiated jets 

Quark- and gluon-initiated jets were selected on a statistical basis based on their sub- 
structure. The integrated jet shape at r = 0.3, ip{r = 0.3), and the subjet multiplicity 
at Hcut = 5 ■ 10~^, nsubjet (2/cut = 5 ■ 10""^), were used to select quark- and gluon-initiated 
jets in the photoproduction and DIS samples. The values r = 0.3 and ?/cut = 5 ■ 10~^ 
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were chosen to be as small as possible to be sensitive to the differences between quark 
and gluons, but large enough to avoid uncertainties due to the detector resolution. The 
different behaviour of these distributions for gluon- and quark-initiated jets is shown in 
Fig. El for samples of Pythia and Herwig generated events. These observables were 
used to classify the jets into: 

• a gluon-enriched sample (broad jets), defined as those jets with ?/j(r = 0.3) < 0.6 
and/or nsubjet(z/cut = 5 ■ 10"^) > 6, and 

• a quark-enriched sample (narrow jets), defined as those jets with ip{r = 0.3) > 0.8 
and/or nsubjct(ycut = 5 ■ 10"^) < 4. 

Non-overlapping ranges were chosen to suppress migration effects. The values for the 
cuts in ip{r = 0.3) and ?^subjct(ycut = 5 ■ 10""^) chosen were a compromise between purity 
and statistics. The purity of the gluon-initiated sample is around 50%, whereas for the 
quark-enriched sample it is around 90%. Table ITUl shows the purities and efficiencies for 
the different MC and selected samples. 

The remaining number of jets after applying the jet-shape, subjet-multiplicity and the 
combination of both selection cuts in the broad and narrow inclusive jet photoproduction 
data sample and in the broad-broad, narrow-narrow and broad-narrow dijet data samples 
selected using the jet-shape method are shown in Table ^2 The same table shows also 
the number of jets in the broad and narrow inclusive jet NC DIS data sample selected 
according to the jet shape. In the next sections, measurements of cross sections are 
presented for samples of jets separated according to their shape and/or subjet multiplicity. 

8.2 Measurements of da/drf^^ 

The differential inclusive-jet cross-section da/drf^^ for photoproduction is shown in Fig.fT^ 
and Table 1121 for samples of broad and narrow jets, separated according to the jet-shape 
selection. The measured cross sections exhibit different behaviours: the rf^^ distribution 
for broad jets increases up to the highest rf^^ value measured, whereas the distribution for 
narrow jets peaks at rf'^^ ~ 0.7. Measurements of da/drf'^^ in photoproduction for samples 
of broad and narrow jets separated according to either the subjet-multiplicity selection 
or a combination of jet shape and subjet multiplicity are shown in Figs. IT^ andlT^ and 
Tables IT^ and ITU respectively. These measurements also exhibit a difference in shape for 
the samples of broad and narrow jets. The same conclusions as in the case of using the 
integrated-jet-shape selection method can be drawn. 

Leading-logarithm parton-shower MC calculations using Pythia, Herwig and Pythia 
MI for resolved plus direct photon processes are compared to the measurements in Figs.fT^ 
andlT^. The same selection method was applied to the jets of hadrons in the MC event 
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samples and the calculations have been normalised to the total measured cross section of 
each sample. The MC predictions provide a good description of the shape of the narrow-jet 
distribution in the data. The shape of the broad-jet distribution in the data is reasonably 
well described by Pythia or Pythia MI, but the prediction of Herwig fails to describe 
this distribution. From the calculation of Pythia (Herwig), the sample of broad jets se- 
lected according to the jet shape is predicted to contain 15(12)% of gg subprocesses in the 
final state and 50(47)% of gq, and a contamination from processes with only quarks in the 
final state of 35(41)%. There is a large contribution from gq final states in the broad-jet 
sample because the partonic cross section for the resolved subprocess q^gp —>■ qg is much 
larger than the cross section for the subprocesses qq — >• gg plus gg gg. The sample of 
narrow jets contains 62(61)% of qq subprocesses and 34(36)% of qg, with a contamination 
of 4(3)% from gg subprocesses. The measured cross section for the broad-jet sample is 
(32.1 ± 0.1)% of the total cross section, whereas the narrow-jet sample is (40.6 ± 0.1)% 
using the jet-shape selection. Pythia (Herwig) predicts 31.5(27.1)% for the broad-jet 
sample and 37.4(44.0)% for the narrow-jet sample. Similar conclusions can be drawn 
from the subjet-multiplicity selection and the combined jet-shape and subjet-multiplicity 
selection. 

Figure ITHk also shows the predictions of Pythia for jets of quarks and gluons sepa- 
rately. These predictions have been obtained without any jet-shape selection and are 
normalised to the data cross sections. The calculation that includes only quark-initiated 
jets gives a good description of the narrow-jet cross section, whereas the calculation for 
gluon-initiated jets provides a reasonable description of the broad-jet cross section. This 
result supports the expectation that the broad (narrow)-jet sample is dominated by gluon 
(quark)-initiated jets. 

Figure and Table ^1 show da/drf'^^ in DIS for samples of broad and narrow jets 
separated according to the jet shape. The two cross sections have the same variation with 
rf'^^, as can be seen from the ratio of the narrow- to the broad-jet cross sections. However, 
the narrow-jet cross section is about five times larger than the broad-jet cross section, 
which shows that the DIS sample is enriched in quark-initiated jets. Since the ratio of the 
cross sections is approximately constant, the quark and gluon content of the final state in 
DIS does not change with rf'^^, as also was concluded from Fig. [7^. The predictions of the 
CDM model are compared to the data in Fig. IT4k and give a good description of the data. 
The predictions of MEPS give a poorer description. Figure shows the same measured 
cross sections together with the calculations of CDM for gluon- and quark-initiated jets; 
no jet-shape selection has been applied in this case. These predictions describe well the 
shapes of the broad- and narrow-jet samples, respectively. 
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8.3 Measurements of da/dE^^^ 



The differential inclusive jet cross-section da/dE-!^', measured in the range 17 < -E;^* < 
95 GeV, is presented in Figs. IT^ and ITSb and Tables El and El for samples of broad 
and narrow jets, separated according to the jet shape, for photoproduction and DIS 
events, respectively. The cross sections for the narrow-jet samples have a harder spectrum 
than that for the broad-jet sample. Figures Efc and ITSH and Tables El and El show 
the da/dE^r^^ cross section for samples of broad and narrow jets separated according 
to the subjet-multiplicity selection and the combined integrated-jet-shape and subjet- 
multiplicity selection in the photoproduction regime. These measurements exhibit the 
same behaviour as in Fig. IT3k . but the cross-over between the broad- and narrow-jet 
distributions takes place at slightly higher i?;^*. The MC calculations using Pythia, 
which have been obtained using the same selection method as for the data, are compared 
to the measurements in Figs. IT^ . c and d. The MC predictions provide a good description 
of the shapes of the data distributions. The predictions of the CDM MC are compared to 
the measurements in Fig.llSb and give a good description of the data. In photoproduction, 
the different E^r^^ spectra exhibited by the narrow- and broad-jet samples are understood 
in terms of the increasing fraction of gluon-initiated jets as -E^* decreases. 



8.4 Measurements of da /d cos 9* 

For samples of broad-broad or narrow-narrow dijet events, only the absolute value of 
cos^* can be determined because the outgoing jets are indistinguishable. The differential 
dijet cross section as a function of | cos^*| has been measured in the range | cos^*| < 0.8 
for dijet invariant masses M^j > 52 GeV for photoproduction. The region of phase 
space in the (M-ij, | cos^*|) plane was chosen in order to minimise the bias introduced 
by selecting jets with E^^^^ > 17 GeV and E^^^'^ > 14 GeV. The measured da/d\ cos^*| 
for the samples of broad-broad dijet events and narrow-narrow dijet events are presented 
in Fig. EK and Table (201 The measured and predicted cross sections were normalised 
to unity at |cos^*| = 0.1. The |cos6'*| distribution for the two samples of dijet events 
increases as | cos 6* \ increases; however they exhibit a different slope. The cross section 
at I cos 6** I = 0.7 for broad-broad dijet events is more than seven times larger than the 
measured value at | cos6'*| = 0.1, whereas for narrow-narrow dijet events, the cross section 
at I cos^*| = 0.7 is only twice as large as at | cos^*| =0.1. 

Calculations using Pythia for broad-broad and narrow-narrow dijet events are compared 
to the data in Fig. IT^ . The predictions from Pythia give a good description of the shape 
of the measured da/d\ cos^*|. Pythia predicts 16% of (^^f-final-state subprocesses, 52% of 
gq and 32% of qq for the broad-broad dijet sample in the kinematic region of this measure- 
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ment. For the narrow-narrow dijet sample, the predictions are: 71% of qq subprocesses, 
28% of qg and 1% of gg. The differences observed in the measured da/d\ cos^*| for the 
two samples are adequately reproduced by the calculations and understood in terms of 
the dominant two-body processes: the resolved subprocess q^gp — ^ qg, mediated by gluon 
exchange for the broad-broad dijet sample and the direct subprocess 7(7 qq, mediated 
by quark exchange for the narrow-narrow dijet sample. 

The sample of photoproduced dijet events with one broad jet and one narrow jet allows 
a measurement of the unfolded da / dcosOl^^^^ cross section. Since in this case the two 
jets can be distinguished, 0"^^^^ refers to the scattering angle measured with respect to 
the broad jet. Figure ITHb and Table E] show the measured dijet cross section as a func- 
tion of cos6'jJj,Q3^j. The measured and predicted cross sections were normalised to unity at 
cos6'^j,Q3^j = 0.1. The dijet angular distribution shows a different behaviour on the negative 
and positive sides; the measured cross section at cosO"^^^^ = 0.7 is approximately twice 
as large as at cos^j^j.^^^^ = —0.7. The calculation from Pythia gives a good description 
of the shape of the measured da / d cos 6^^^^^. The predictions of Pythia for the partonic 
content are: 52% of qg subprocesses, 4% of gg and 44% of qq. The observed asymmetry 
is adequately reproduced by the calculation and is understood in terms of the dominant 
resolved subprocess q^gp — * qg. The cosOlj.^^^ distribution for this subprocess is asym- 
metric due to the different dominant diagrams in the regions cos^j^j.^^^^ — ^ ±1: t-channel 
gluon exchange {cosO^^^^^ ^ +1) and u-channel quark exchange {cosO^^^^^ — > — 1). 



8.5 Measurements of da/dM^^ and da/dx°^^ 

The photoproduction differential dijet cross section as a function of M^j has been mea- 
sured in the range 52 < M^j < 123 GeV for |cos6'*| < 0.8. The measured da/dM^^ 
for the samples of broad-broad and narrow-narrow dijet events are presented in Fig. ITHh 
and Table 1221 The measured da/dM^^ cross sections decrease as MJj increases, but the 
distribution for the narrow-narrow dijet sample exhibits a harder spectrum, as was also 
seen for the inclusive jet cross section as a function of i?;^*. The MC calculations from 
Pythia are compared to the data in Fig. IT^ and give a good description of the shape 
of the measured da/dM^K The different shape in both cross sections is understood in 
terms of the dominant two-body processes: the broad-broad dijet sample is dominated by 
the resolved subprocess q-ygp qg and the narrow-narrow dijet sample is dominated by 
the direct subprocess 'jg — > qq. Direct processes reach larger values of M-i-i than those of 
resolved since the full incoming-photon energy is available at the hard interaction. 

Resolved and direct processes can be separated by using the x'^^ variable, which is defined 
as 
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Resolved and direct processes populate different regions in x?l°'^, with the direct processes 
concentrated at high values. The dijet cross section as a function of x"^^ is presented 
in Fig. ITBH and Table ESI and is reasonably well described by the MC predictions of 
Pythia. The cross section for the broad-broad dijet sample is approximately constant 
as a function of x'^^ whereas the cross section for the narrow-narrow dijet sample peaks 
at high values. The shape of the distribution for the broad-broad (narrow-narrow) dijet 
events is consistent with the dominance of resolved (direct) processes. 



9 Determination of as 

The measured {i/j{r = 0.5)) for E^r^^ > 21 GeV in DIS was used to determine as{Mz) 
using a method similar to one presented previously [27]. The NLO QCD calculations 
were performed using the program DiSENT with three different MRST99 sets of proton 
PDFs, central, MRST99ii and MRST99TT; the value of ^.(M^) used in each partonic 
cross-section calculation was that associated with the corresponding set of PDFs. The 
as{Mz) dependence of the predicted {ilj{r = 0.5)) in each bin i of E^r^^ was parameterised 
according to 

mr = 0.5))(a,(Mz))], = Cl + Cla^iMz), 

where C{ and C2 were determined from a fit by using the NLO QCD calculations 
corrected for hadronisation and Z'^-exchange effects. Finally, a value of as{Mz) was 
determined in each E^^^ region as well as from all the data points by a fit. 

The values of as{Mz) as determined from the measured {ip{r = 0.5)) in each region 
of Ej^^ are shown in Fig. El and Table |^ Taking into account only the statistical 
uncertainties, the value of as{Mz) obtained by combining all the i?;^* regions is as{Mz) = 
0.1176 ±0.0009 (stat.). 

The uncertainties on the extracted value of as{Mz) due to the experimental systematic 
uncertainties were evaluated by repeating the analysis above for each systematic check 
described in Section 16.31 The total experimental systematic uncertainty on the value of 
as{Mz) is ^as{Mz) / Ois^Mz) = -2^%- The main contribution to the positive (negative) 
systematic uncertainty comes from the uncertainty in the jet energy scale (scattered- 
electron identification). 

The following sources of theoretical uncertainties on the extracted value of as{Mz) were 
considered: 

• terms beyond NLO were estimated by varying fiji between Q/2 and 2Q and keeping 
/ii? fixed at Q; this results in a variation of Aas(M^) = lo;oo7o; 
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• the uncertainty on the modeUing of the parton shower was estimated by using the 
MEPS model to calculate the parton-to-hadron correction factors; this results in a 
variation of Aa,(M^) = 0.0018; 

• the uncertainty in the choice of was estimated by using /i^j = E^r^^ instead of Q and 
fip was set to Q; this results in a variation of Aas{Mz) = 0.0003; 

• the uncertainty in the NLO QCD calculations due to the uncertainties in the pro- 
ton PDFs was estimated by repeating the calculations using 40 additional sets from 
CTEQ6 [53]; this results in a variation of Aa,(Mz) = 0.0002; 

• the uncertainty of the calculations in the value of fip was estimated by repeating the 
calculations with fip = Q/2 and 2Q; this results in a variation of Aas{Mz) = 0.0001. 

These uncertainties were added in quadrature and give a total theoretical uncertainty 
of Aas{Mz)/a,{Mz) = As a cross-check of the extracted value of as{Mz), the 

fit procedure was repeated by using the five sets of the CTEQ4 "A-series" , resulting in 
as{Mz) = 0.1178 =b 0.0009, in very good agreement with the central value determined 
above. As a consistency check, the whole procedure was repeated for {iplr = 0.4)) and 
(z/;(r = 0.6)), giving values of a,(Mz) = 0.1158 ± 0.0008 and a,(Mz) = 0.1193 ± 0.0010, 
respectively, which are compatible with the value determined from (?/'(r = 0.5)). The 
determination of as{Mz) was also repeated using the calculations from the Disaster++ 
program; this gave as{Mz) = 0.1166 ± 0.0009, which is compatible with the value quoted 
above. 

The value of as{Mz) as determined from the measured {ip{r = 0.5)) is therefore 

as{Mz) = 0.1176 ±0.0009 (stat.) t'.^l (exp.) 1°;°°?^ (th.). 

This result is in agreement with recent determinations using measurements of jet produc- 
tion in DIS [1,4,6,8,27] and photoproduction [5] and with the current world average of 
0.1183 ± 0.0027 [54]. This determination of as has experimental uncertainties as small as 
those based on previous measurements. However, the theoretical uncertainty is large and 
dominated by terms beyond NLO. Further theoretical work on higher-order contributions 
would allow an improved determination of from the integrated jet shape in DIS. 

10 Summary and conclusions 

Measurements of the mean integrated jet shape and mean subjet multiplicities for in- 
clusive jet photoproduction and DIS at a centre-of-mass energy of 318 GeV using the 
data collected by ZEUS in 1998 to 2000, which correspond to an integrated luminosity 
of 82.2 pb~^, have been presented. The measurements refer to jets identified with the kx 
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cluster algorithm in the longitudinally invariant inclusive mode in the laboratory frame 
and selected according to E^r^^ > 17 GeV and —1 < rj-'^^^ < 2.5. The measurements are 
given in the kinematic region defined by < 1 GeV^ and 142 < W^p < 293 GeV for 
photopro duct ion and > 125 GeV^ for DIS. The jet shape broadens (narrows) and the 
mean subjet multiplicity increases (decreases) as rj^^^ (E^t^) increases in photoproduction. 
The observed broadening of the jet shape and the increase of the mean subjet multiplicity 
as rf'^^ increases are consistent with an increase of the fraction of gluon-initiated jets. In 
DIS, the data show no significant dependence with rf^^ and a moderate dependence with 
E^r^^. Leading-logarithm parton-shower MC models for photoproduction and NLO QCD 
calculations for DIS give a good description of the data. 

Measurements of differential inclusive jet and dijet cross sections in the photoproduc- 
tion and DIS regimes separated into broad and narrow jets according to their internal 
structure have been presented. Leading-logarithm parton-shower MC models give a good 
description of the data. The inclusive jet cross-sections da/dr]^^^ and da/dE^r^^ for broad- 
and narrow-jet samples show the expected behaviour for samples enriched in gluon- and 
quark-initiated jets, respectively. The dijet cross section as a function of | cos 6^*1, mea- 
sured in the range | cos^*| < 0.8 and integrated over M^j > 52 GeV, displays for broad- 
broad dijets a behaviour consistent with that expected for a sample enriched in processes 
mediated by gluon exchange. Narrow-narrow dijets, however, show a behaviour consis- 
tent with a sample enriched in processes mediated by quark exchange. The dijet cross- 
section da / d cos 9^j.^^^, measured in the region — 0.8 < cos 6*^;^,^^^^ < 0.8 and integrated over 
> 52 GeV, for a sample of events with broad-narrow dijets, exhibits a large asymme- 
try consistent with the expected dominance of gluon (quark) exchange as cos9l^^^^ — > +1 
{cos 91^^^ — > — 1). The dijet cross section as a function of M^j and x'^^ for the sample 
of broad-broad dijets shows a behaviour consistent with the dominance of the resolved 
Q-ydp subprocess, whereas the sample with narrow-narrow dijets is consistent with 

the dominance of the direct subprocess 7(7 qq. 

The measurements of the mean integrated jet shape in DIS have been used to extract 
a value of as{Mz) by comparing to the predictions of NLO QCD as a function of E^^^ . 
The calculations reproduce the measured observables well, demonstrating the validity 
of the description of the internal structure of jets by pQCD. The value of as{Mz) as 
determined by fitting the NLO QCD calculations to the measured mean integrated jet 
shape {i){r = 0.5)) for E^?* > 21 GeV is 

as{Mz) = 0.1176 ±0.0009 (stat.) 1°;°°°^ ^^^^^^ +0.0091 ^^^ ^ 
This value is in good agreement with the current world average. 
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Table 1: Measured mean integrated jet shape corrected to the hadron level for 



jets in photoproduction with E';^* 
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Table 2: Measured mean integrated jet shape corrected to the hadron level for jets 
in photoproduction with — 1 < rj^^^ < 2.5 in different E^^^ regions. The statistical 
and systematic uncertainties are also indicated. 
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Table 3: Measured mean integrated jet shape corrected to the hadron level for jets 
in photoproduction with —1 < rj^^^ < 2.5 in different E^^^ regions. The statistical 
and systematic uncertainties are also indicated. 
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Table 4: Measured mean integrated jet shape corrected to the hadron level and 
for electroweak radiative effects for jets in DIS with E^f:^ > 17 GeV in different rj^^^ 
regions. The statistical and systematic uncertainties are also indicated. 
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Table 5: Measured mean integrated jet shape corrected to the hadron level and 
for electroweak radiative effects for jets in DIS with —1 < rf'^^ < 2.5 in different 
E'^* regions. The statistical and systematic uncertainties are also indicated. 
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Table 6: Measured mean integrated jet shape corrected to the hadron level and 
for electroweak radiative effects for jets in DIS with —1 < rf'^^ < 2.5 in different 
regions. The statistical and systematic uncertainties are also indicated. 
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Table 7: Measured mean subjet multiplicity corrected to the hadron level for jets 
in photoproduction with E^f^ > 17 GeV in different 77-''^* regions. The statistical and 
systematic uncertainties are also indicated. 
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Table 8: Measured mean subjet multiplicity corrected to the hadron level for jets 
in photoproduction with —1 < r]^'^^ < 2.5 in different E^rp" regions. The statistical 
and systematic uncertainties are also indicated. 



33 



Vcut 



0.0005 
0.001 
0.003 
0.005 
0.01 
0.03 
0.05 
0.1 



Vcut 



0.0005 
0.001 
0.003 
0.005 
0.01 
0.03 
0.05 
0.1 



('^subjet (ycut)) ± Stat. ± syst. 



35 < ^J?* < 41 GeV 



4.152 ±0. 
3.264 ±0. 
2.272 ±0. 
1.885 ±0. 
1.514 ±0. 
1.226 ±0. 
1.143 ±0, 
1.052 ±0. 



-0.030 
-0.090 
-1-0.022 
-0.076 
-1-0.017 
-0.063 
-1-0.015 
-0.025 

m 9-1-0.028 
'-'-'-^-0.012 

QQq-HO.014 
"-"-"^-O-Oll 
-1-0.007 
-0.015 



027: 
022 
016 
014 



007 



nn e:-H0.005 
"-'"-'^-0.005 



41 < ^j?* < 47 GeV 



1.805 ±0.022 



4.000 ± 0.043+1^:^35 

3.128 ± 0.0351^23 
2.148 ± 0.026l|]:[|^^ 

-1-0.023 
-0.036 

1.484 ± 0.020l[]:[|^^ 
1.202 ± 0.015l°:°^2 

1.113 ± 0.01 il|]:[J}^ 

1.045 ± 0.007l|]:[}?^ 



('^subjct(ycut)) ± Stat. ± syst. 



47 < ^j?* < 55 GeV 



3.702 ± 
2.872 ± 
2.006 ±0 
1.725 ±0 
1.411 ±0 
1.186 ±0 
1.134 ±0 
1.028 ±0 



Qf^o-HO.063 
"-'"'-'-0.088 
nc:i+0.051 
"-'^-'--0.052 
nQ7-l-0.037 
'-''J ' -0.038 
r|qq-l-0.035 
'-''J'J-0.034 
Q^q-^0.030 
"-'■^^-0.030 
r|9q-l-0.023 
'-'^'J-0.033 

ni O-I-0.018 
U-l-O-0.025 

QQq-HO.012 



55 < E}^' 



< 71 GeV 



3.621 ±0 
2.726 ± 
1.971 ±0 
1.722 ±0 
1.459 ± 
1.268 ±0 
1.172 ±0 
1.081 ±0 



Qq4+o.095 

"-'^^-0.302 

n7i -H0.073 
-'--0.074 

np;p;-H0.055 
UJJ-0.082 

nc;Q-H0.054 
'-'^'-'-0.079 

'-'^'-'-0.068 
nqo-l-0.038 
"-''-'°-0.060 
nqi -1-0.031 
'-''-'-'--0.072 
m q-f 0.020 
'-'-'-^-0.061 



Table 9: Measured mean subjet multiplicity corrected to the hadron level for jets 
in photoproduction with —1 < r]^^^ < 2.5 in different E^rp" regions. The statistical 
and systematic uncertainties are also indicated. 
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Pythia 


Herwig 


CDM 


MEPS 


jet shape (%) 


gluons 
quarks 


57/48 
84/51 


50/43 
85/55 


43/47 
93/65 


33/65 
97/57 


subjet multiplicity (%) 


gluons 
quarks 


56/63 
90/23 


53/56 
92/32 




combination (%) 


gluons 
quarks 


61/43 
90/22 


55/37 
93/30 





Table 10: Purity /efficiency of gluon-initiated jets in the broad-jet sample and 
of quark-initiated jets in the narrow-jet sample from the Pythia and Herwig 
photoproduction Monte Carlo generated samples and from the CDM and MEPS 
NC DIS Monte Carlo generated samples. 



jet shape 





photoproduction 


NC DIS 


broad 
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narrow 
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broad-broad 
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broad-narrow 
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subjet multiplicity 





photoproduction 


broad 


55750 


narrow 


54654 


combination 




photoproduction 


broad 


37213 


narrow 


43928 



Table 11: Number of jets (events) in the inclusive-jet (dijet) photoproduction 
and NC DIS samples selected according to various methods. 
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( .i HZ U.O_2.6-0.9 

60.3 ±1.11}^:^!^:^ 

130.0 ±1.6™^:^ 

240.1 ±2.5^t6-n:6 
259.1 ih 2.5^]^]^ y]'^]^2.2 
270.8 ± 2.5^}^:?^}2.7 


20.9 ± 0.7+};^!};^ 

137.8 ±1.7™|5^ 

07 A Q -1- p:+26. 5+10.8 
Z/^.O in 4-11.1 

349.0 ± 2.8li:^l}J;? 
299.3 ± 2.7+^;°l|^ 

227.9 ± 2.2tlXt^S 
151.7 ±1.8lt}l5:4 



Table 12: Measured differential ep cross-section da/drj^'^^ for inclusive jet photo- 
production with E^rf^ > 17 GeV. The jets have been selected according to their shape 
into broad and narrow jets. The statistical and systematic uncertainties —not asso- 
ciated with the absolute energy scale of the jets— are also indicated. The systematic 
uncertainties associated to the absolute energy scale of the jets are quoted separately. 





da/drf^^±. stat. ± syst. ± -E^r* — scale (pb) 


broad jets 


narrow jets 


-1.0- -0.5 
-0.5-0.0 
0.0-0.5 
0.5-1.0 
1.0-1.5 
1.5-2.0 
2.0 - 2.5 


8.1±0.3t|ni 

79 1 _|_ 1 1 +14.8+4.4 
( z.i HZ -'-•-i-_i4 8_5.3 

153.1 ± 1-6123.6-7.9 

oi o 7 _|_ 9 O+13.0+9.6 
ZiO. t HZ ^•0_i3 0-10.8 

261.4 ±2.8^1™ 

300.7 ±2.9l22;3+12.7 

315.5 ± 2.6li:{llf, 


8.9 ± 0.3l|^l°| 
63.8 ±0.9126-5+3.1 

132.8 ±1.4lti^^:^ 

171.9 ±1.611^^1^:° 

137.4 ±1.2li:ilt3 
62.4±0.8l™.2 



Table 13: Measured differential ep cross-section da/drj^^^ for inclusive jet photo- 
production with E^rf^ > 17 GeV. The jets have been selected according to their subjet 
multiplicity into broad and narrow jets. The statistical and systematic uncertain- 
ties — not associated with the absolute energy scale of the jets— are also indicated. 
The systematic uncertainties associated to the absolute energy scale of the jets are 
quoted separately. 
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da/dr]^^^± stat. ± syst. ± i?;^* — scale (pb) 


broad jets 


narrow jets 


-1.0 - -0.5 
-0.5 - 0.0 
0.0-0.5 
0.5 - 1.0 
1.0 - 1.5 
1.5-2.0 
2.0-2.5 


^. / in v.o_i j_Q Q 

44.0 ± 0.9tTi-l' 
96.7 ±1.4116.4+5.6 

139.9 ii.g^lts-?:' 

182.0 ± 2.3tlitli 
213.5 ±2.4l-^t^°- 
240.5 ± 2.4l^0;3+n-3 


8.5 ± 0.3+1^1°:^ 

58.8 ± 0.9l^2;}+2;8 

121.7 ±lAtir,tti 

160A±1.6tir,ni 
126.5 ±1.3t||t|? 
94.2 ±1.1^11^+2.9 

55.2 ± O.Stltitli 



Table 14: Measured differential ep cross-section da/drj^'^^ for inclusive jet pho- 
toproduction with i?^* > 17 GeV. The jets have been selected according to the 
combination of the jet shape and subjet multiplicity into broad and narrow jets. 
The statistical and systematic uncertainties —not associated with the absolute en- 
ergy scale of the jets— are also indicated. The systematic uncertainties associated 
to the absolute energy scale of the jets are quoted separately. 





da/rfr/j*^* ± stat. ± syst. ± -E;^* — scale (pb) 


broad jets 


narrow jets 


-1.0 - -0.5 
-0.5 - 0.0 
0.0-0.5 
0.5 - 1.0 
1.0 - 1.5 
1.5-2.0 
2.0-2.5 


« _L Q 4+0.9+0.4 

21.1±0.7+°ill-° 
41.4±1.0li^ll| 

51.4 ±1.11^:^111 

57.3 ± 1.311™ 
57.5±l.ll|riS 

49.4 ± 1.01^:111:^ 


26.4 ± 0.916:111:5^ 

95.4 ±1.51^:^11^ 

181.6 ±2.iil^:«ll:^ 

236.1 ±2.511^^11:6^ 
252.3 ± 2.7l|iltl 
256.3 ± 2.61^:11^^ 

211.2 ±2.4lt°ltii 



Table 15: Measured differential ep cross-section da/dtf^^ for inclusive jet pro- 
duction in DIS with E^^^ > 17 GeV. The jets have been selected according to their 
shape into broad and narrow jets. The statistical and systematic uncertainties —not 
associated with the absolute energy scale of the jets— are also indicated. The sys- 
tematic uncertainties associated to the absolute energy scale of the jets are quoted 
separately. 
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E}^' (GeV) 


da/dE^r^^± stat. ± syst. ± E'^* — scale (pb) 


broad jets 


narrow jets 


17. - 21. 
21. - 25. 
25. - 29. 
29. - 35. 
35.-41. 
41. -47. 
47. - 55. 
55. - 71. 
71. - 95. 


102.4 ±0.5l|™ 
27.2 ± 0.3l}il}-^ 
8.78±0.16l°i^l°:i 
3.02 ± 0.08l°:|;°;i 

0.774 ±0.040l°:°^^ini^ 

U.ZOi in U.UZZ_Q Q27_o,020 

0.064 ±0.0101°:°}^::^ 

n nnnc;n+o-ooo88+o.ooo98+o.oooo6 

U.UUUOU„o.00032-0.00001 -0.00009 


102.3 ± 0.5tlitli 
42.7±0.3;|};j;^ 
18.9 ± 0.2t\it',i 

7 oo _|_ n 1 9+0.43+0.41 
' -o^ =■= '-'•-'-^-0.42-0.33 

2.93 ± 0.07l°-}^+°-J3 

1.094 ± OMGtrollt'ofsl 
0.447 ±0.026l°:°-l°;°i 
0.068 ± 0.006+°:°^°!°;:^ 

0.0102 ±0.00271°X1™ 



Table 16: Measured differential ep cross-section da /dE^^ for inclusive jet pho- 
toproduction with —1 < rj^^^ < 2.5. The jets have been selected according to their 
shape into broad and narrow jets. The statistical and systematic uncertainties —not 
associated with the absolute energy scale of the jets— are also indicated. The sys- 
tematic uncertainties associated to the absolute energy scale of the jets are quoted 
separately. 



E}^' (GeV) 


da/dE^rf^± stat. ± syst. ± E'y* — scale (pb) 


broad jets 


narrow jets 


17. - 21. 
21. - 25. 
25. - 29. 
29. - 35. 
35.-41. 
41.-47. 
47. - 55. 
55. - 71. 
71.-95. 


19.6 ±0.2^1:61°] 

7 7« _L n 1 c:+0.75+0.31 

o cr7,Q 1Q+0.30+0.16 
O.OI in U.iU_o.3o_o.l3 

1.70 ±0.061°:°™ 
0.714 ±0.0381°:°^^!°;°^^ 
312 ± n26+°-°^i+°-°°^ 

0.151 ±o.oi5i°;°}^;°:°Si 

0.0349 ± 0.00521°:°°^^!°:°°^^ 
0.0099 ± 0.0024l°:°°J^l°:°°- 


67.3 ± 0.5liJ+|;0 
35.2 ± 0.3liJ10:? 
20.7±0.3+i-}l°t 

10.81 ±o.i5irS:l 
5.36±o.iii°:^^l°:i? 
3.02 ± o.o8l°:i^l°:°? 

1 468 ± n4R+o.035+o.073 

i.^UO m U.U^U_o.o35-0.052 
S72 ± 021+0-035+0.017 

llfic^ -U QQ74+0.0081+0.0089 

U.iiOO HI U.UU r^„Q QQg;^_Q QQ55 



Table 17: Measured differential ep cross-section da/dEi^ for inclusive jet pro- 
duction in DIS with —1 < 77-''^* < 2.5. The jets have been selected according to their 
shape into broad and narrow jets. The statistical and systematic uncertainties —not 
associated with the absolute energy scale of the jets— are also indicated. The sys- 
tematic uncertainties associated to the absolute energy scale of the jets are quoted 
separately. 
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E'r^' (GeV) 


da/dE^r^^± stat. ± syst. ± -E^* — scale (pb) 


broad jets 


narrow jets 


17. - 21. 
21. - 25. 
25. - 29. 
29. - 35. 
35.-41. 
41.-47. 
47. - 55. 
55. - 71. 
71.-95. 


116.6 ±0.6^1i+5;i 

32.7±0.3l°-^l}-^ 
10.91 ±0.181™J^ 

3.33 ± omt'oT-'oil 

0.868 ± 0.0401°:°?^+°:°^° 
n 997 -u n noi +0.030+0. oie 

V.ZZI IC U.UZi_o.o30-0.018 

0.084 ±0.0101°:°!^^:°°^ 

0.0133 ±0.0028l°;°^il°:rr 
n nnm fi+o. 00028+0.00017+0. 00000 
u.uuuiu_gQooio_o.oooo4-o. 00002 


44.7 ±0.3+}^!+};^ 
19.7 ±0.2^^6-0:? 

9 + 1+2.4+0.3 
J.Z HZ U.l„2.4-0.5 

4 38 + 08+0-65+0.22 

^.00 in U.UO_Q g4_Q 20 

1.74±0.05l°-n°^ 

0.70 ± omtrilTol 

n Q97 -u n noT +0.015+0.028 

V.OZI in U.UZ1_Q Q]^5„Q Q2g 

0.055 ± Om5t'o'ollt'oZ 
0.0104 ±0.0028^^°!^°:^^ 



Table 18: Measured differential ep cross-section da /dE^^ for inclusive jet pho- 
toproduction with —1 < rf^^ < 2.5. The jets have been selected according to their 
subjet multiplicity into broad and narrow jets. The statistical and systematic un- 
certainties —not associated with the absolute energy scale of the jets— are also 
indicated. The systematic uncertainties associated to the absolute energy scale of 
the jets are quoted separately. 





(GeV) 


17. 


- 21. 


21. 


-25. 


25. 


- 29. 


29. 


-35. 


35. 


-41. 


41. 


-47. 


47. 


-55. 


55. 


- 71. 


71. 


-95. 



da/dEj:^± stat. ± syst. ± E'y* — scale (pb) 



broad jets 



narrow jets 



82.7±0.5t};Jlt2 
21 4 ± 2+°-^+^-2 
6 72 + 14+0.14+0.36 



0.10+0.16 
0.10-0.17 
-0.062+0.023 
-0.063-0.029 



2.18 ±0.07; 
0.505 ±0.032 
0.138 ±0.017;°:°1}1°;°}J 
0.044 ± 0.0081°:°°^!°;°°^ 

0.0071 ± o.oo2il°:°°-^n°J° 



40.5±0.3tJ}:Jl|i 

18.2 ± 0.2tt'otU 
8.7±0.ll|2+o.2 

4.12 ±0.081°:^™ 

i.66±o.o5l°:rn^ 



0.67 ±0.03 



+0.12+0.02 

0.12-0.05 



Q 017 I Q Q91+0.013+0.030 
U.Oi/ in U.UZi_Q Q^3_Q Q26 

0.051 ±o.oo5i°:°^n:S°^ 
0.0091 ±o.oo27i°::-^n^SS 



Table 19: Measured differential ep cross-section da /dE^^^ for inclusive jet pho- 
toproduction with —1 < 77-^'^* < 2.5. The jets have been selected according to the 
combination of the jet shape and subjet multiplicity into broad and narrow jets. 
The statistical and systematic uncertainties —not associated with the absolute en- 
ergy scale of the jets— are also indicated. The systematic uncertainties associated 
to the absolute energy scale of the jets are quoted separately. 
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COS 9* 


da / d\cos6*\ ±. stat. ± syst. ± i?;^* — scale (pb) 


broad-broad dijets 


narrow- narrow dijets 


0.0-0.2 

0.2 - 0.4 
0.4-0.6 
0.6 - 0.8 


3.93 ± 0.49l°:i^°:i 

7 Q _|_ n 7+1.4+0.5 
/ .O in U. / _i.4_o.4 

10.4±0.8l};il?:^ 
29.9 ± 1.2+1^+?:° 


43.9±1.6li™ 

49.0± 1.7lt2t2:? 
67.9 ± 2.Qt\ltll 
86.3 ± 2^t^l 



Table 20: Measured differential ep cross-section da /d\ cos 9*\ for dijet photo- 
production with iVfjj > 52 GeV. The events have been selected according to their 
shape into broad-broad and narrow-narrow dijets. The statistical and systematic 
uncertainties —not associated with the absolute energy scale of the jets— are also 
indicated. The systematic uncertainties associated to the absolute energy scale of 
the jets are quoted separately. 





da / d cos OIj.^^^ ± stat. ± syst. ± .Ey* — scale (pb) 


cos ^broad 


broad-narrow dijets 


-0.8 - -0.6 


30.3±1.2lt2^}1 


-0.6 - -0.4 


16.5 ± 0.9tlit'oi 


-0.4 - -0.2 


14.5 ± 0.9t|°in 


-0.2-0.0 


10.3 ± O.Ttlit'o'r 


0.0-0.2 


12.8 ± 0.8l}:^^|]j 


0.2 - 0.4 


i7.o±i.ol|ll?:S 


0.4-0.6 


31.7±1.4l}-^l}-^ 


0.6-0.8 





Table 21: Measured differential ep cross-section da /d cos 6^^.^^^ for dijet photo- 
production with iVfjj > 52 GeV. The events have been selected according to their 
shape into broad-narrow dijets. The statistical and systematic uncertainties —not 
associated with the absolute energy scale of the jets— are also indicated. The sys- 
tematic uncertainties associated to the absolute energy scale of the jets are quoted 
separately. 
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MJj (GeV) 


da/dAP^± stat. ± syst. ± iJ^* — scale (pb) 


broad-broad dijets 


narrow-narrow dijets 


52. - 57. 

57. - 65. 

65.- 83. 
83. - 103. 
103. - 123. 


1 006 ± 047+0-038+0.060 
i.UUU HZ U.U^r „Q Q3g„Q 

n Qfi7 -U n 099+0-035+0.030 
V.OOl in U.UZZ_Q Q35_Q 

0.105 ±0.008tHrH^^ 

0.0109 ± 0.0023;°:°°}^^°:°°}? 

n nni k;+o.ooi2+o.ooi6+o.ooo3 

U.UUiO_o.ooiO-0.0003-0.0001 


3.06 ± 0.08t'oiV-ril 

1.842 ±0.0531°:°!;°:°^? 

0.749 ± 0.022lH?^tH^? 

0.210 ±o.oiit°:°°^tnS^ 
0.0574 ±o.oo59l°::|l°::?^ 



Table 22: Measured differential ep cross-section da/dM^^ for dijet photoproduc- 
tion with \ cos^*| < 0.8. The events have been selected according to their shape into 
broad-broad and narrow-narrow dijets. The statistical and systematic uncertain- 
ties — not associated with the absolute energy scale of the jets— are also indicated. 
The systematic uncertainties associated to the absolute energy scale of the jets are 
quoted separately. 



™obs 


da/dx'^^± stat. ± syst. ± E"^* — scale (pb) 


broad-broad dijets 


narrow-narrow dijets 


0.0-0.2 
0.2 - 0.4 
0.4-0.6 
0.6-0.8 
0.8-1.0 


31.6±1.5tl-ltl-l 
79.2 ± 2.3ttlttl 
66.0 ± 2.0t\-ltli 
71.9±2.0l}}:^l|^ 

70 n -h 1 q+20.6+7.0 
' ^-"J =■= -'-•^-20.3-8.9 


QQ + Q 0+0.9+0.3 
y.u in u.o„o.8-o.3 

48.7± 1.7^^:^^?j 
78.1 ±2.21^3+0.8 

131.1 ±2.8l}|^l?i 
459.6 ± 



Table 23: Measured differential ep cross-section da /dx'^^^ for dijet photoproduc- 
tion. The events have been selected according to their shape into broad-broad and 
narrow-narrow dijets. The statistical and systematic uncertainties —not associated 
with the absolute energy scale of the jets— are also indicated. The systematic un- 
certainties associated to the absolute energy scale of the jets are quoted separately. 
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(GeV) 


as{Mz) ± stat. ± syst. ± th. 


22.9 
26.8 
31.6 
37.7 
43.7 
50.5 
58.5 
79.3 


n 1 1 77 + n nni q+o.ooo5+o.oo96 

V.LII 1 IC U.UUiO_o.o027-0.0076 

0.1177 ±o.ooi8l°:r^nS?^ 

0.1179 ±0.002linS-l°:rr? 
0.1188 ±0.0032l°:0°^°l™ 

0.1134 ±0.0045l°Sl°::n 

0.1165 ±0.0058l°:r}^S 
0.1126 ±0.0068l°:°Jil°:ll 

0.120 ±0.0141°:°?^!°:°°^ 



Table 24: The as{Mz) values determined from the QCD fit of the measured 
{'ijj{r = 0.5)) as a function of E^r^^ in DIS. The statistical, systematic and theoretical 
uncertainties are also indicated. 
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Figure 1: Measured mean integrated jet shape corrected to the hadron level (dots), 
{ip{r)) , for jets in photoproduction with E^^^ > 17 GeV in different rj^^^ regions. The 
error bars, which are typically smaller than the dots, show the statistical and system- 
atic uncertainties added in quadrature. For comparison, the predictions o/Pythia 
including resolved plus direct processes for quark (dot-dashed lines), gluon (dashed 
lines) and all (solid lines) jets are shown. The open circles show the fractional 
difference of the data to the predictions of Pythia for all jets. 
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Figure 2: Measured mean integrated jet shape corrected to the hadron level (dots), 
{tp{r)), for jets in photoproduction with E}^^ > 17 GeV in different rj"^^^ regions. For 
comparison, the predictions o/Herwig (dashed lines) and Pythia MI (dot-dashed 
lines) including resolved plus direct processes are shown. The open circles show the 
fractional difference of the data to the predictions o/Herwig. Other details are as 
in the caption to Fig. QJ 
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Figure 3: Measured mean integrated jet shape corrected to the hadron level (dots), 
{ip{r)), for jets in photoproduction in the range —1 < 77-'^* < 2.5 in different E^^^ 
regions. For comparison, the predictions of Pythia including resolved (dashed 
lines), direct (dot-dashed lines) and resolved plus direct processes (solid lines) are 
shown. Other details are as in the caption to Fig. 
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Figure 4: Measured mean integrated jet shape in photoproduction corrected to the 
hadron level at a fixed value of r = 0.5 (dots), {il){r = 0.5)), as a function of (a) 
ri^°^ with > 17 GeV and (h) E^r^^ with —1 < 77-^'^* < 2.5. Other details are as in 
the captions to Figs. Q and\^ 
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Figure 5: Measured mean integrated jet shape corrected to the hadron level and for 
electroweak radiative effects (squares), {il}{r)) , for jets in DIS with E^^^ > 17 GeV in 
different 77-'°* regions. For comparison, NLO predictions corrected for hadronisation 
and -exchange effects (solid lines) are shown. Other details are as in the caption 
to Fig. Ql 
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Figure 6: Measured mean integrated jet shape corrected to the hadron level and 
for electroweak radiative effects (squares), {ip{r)), for jets in DIS in the range 
— 1 < 77-'*^* < 2.5 in different E^r^^ regions. Other details are as in the captions to 
Figs. U\ and{^ 
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Figure 7: Measured mean integrated jet shape in DIS corrected to the hadron 
level and for electroweak radiative effects at a fixed value of r = 0.5 (squares), 
{i/j{r = 0.5)), as a function of (a) rj^^^ with E^r^^ > 17 GeV and (h) E^^'^ with 
— 1 < r/j'^* < 2.5. Other details are as in the caption to Fig. [3 
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Figure 8: Measured mean integrated jet shape corrected to the hadron level at a 
fixed value of r = 0.5 for DIS (squares) and photoproduction (dots), {ip{r = 0.5)), 
as a function of (a) rj^^^ with E^^^ > 11 GeV and (h) E^r^^ with —1 < 77^'^* < 2.5. 
The predictions for gluon-initiated (lower hatched areas) and quark-initiated (upper 
hatched areas) jets are also shown. The bounds of each hatched area are given by 
the predictions of CDM and Pythia. Other details are as in the caption to Fig. O 
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Figure 9: Measured mean suhjet multiplicity corrected to the hadron level (dots), 
('^subjetd/cut)), for jets in photoproduction with E^^^ > 17 GeV in different r]^°^ 
regions. Other details are as in the caption to Fig. {J\ 
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Figure 10: Measured mean subjet multiplicity corrected to the hadron level (dots), 
('^subjetd/cut)); for jcts in photoproduction in the range —1 < rj^^^ < 2.5 in different 
E^r^^ regions. Other details are as in the caption to Fig. 
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Figure 11: Measured mean subjet multiplicity in photoproduction corrected to 
the hadron level at a fixed value of ycnt = 10~^ (dots), (nsubjct(ycut = 10"^), as a 
function of (a) V'^* with E^^^ > 17 GeV and (h) E}^^ with -1 < < 2.5. Other 
details are as in the caption to Fig. ^ 
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Figure 12: (a) The predicted integrated jet shape distribution at r = 0.3 and (c) 
the predicted subjet multiplicity distribution at ?/cut = 5 ■ 10~^ at the hadron level 
for samples of gluon- ( shaded histograms ) and quark-initiated (hatched histograms ) 
jets simulated using the program Pythia; (b) and (d) show the same distributions 
for samples of Herwig. 
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Figure 13: Measured differential ep cross-section da/drj^'^^ for inclusive jet pho- 
toproduction with E^r^^ > 17 GeV in the kinematic region defined by < 1 GeV^ 
and 142 < W.yp < 293 GeV. The jets have been selected according to (a,b) their 
shape, (c) subjet multiplicity or (d) a combination of both in broad jets (dots) and 
narrow jets (open circles). The thick error bars (not visible) represent the statistical 
uncertainties of the data, and the thin error bars show the statistical and systematic 
uncertainties —not associated with the uncertainty in the absolute energy scale of 
the jets, shown as a shaded band— added in quadrature. The calculations o/Pythia 
for resolved plus direct processes separated according to the same criteria as in the 
data are included in (a,c,d) (solid lines). In (a), the calculations 0/ Pythia for 
gluon (dashed line) and quark (dot-dashed line) jets are also included. In (b), the 
calculations of Herwig (dashed lines) and Pythia MI (dot-dashed lines) with 
the same selection as in the data are included. The MG calculations have been 
normalised to the total measured cross section of each type. In (c) and (d), the 
measurements and predictions for narrow jets have been plotted at 7]^°^ + 0.05 for 
clarity of presentation. 
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Figure 14: Measured differential ep cross-section da/drf'^^ for inclusive jets in 
DIS with E'y* > 17 GeV in the kinematic region defined by > 125 GeV^. The 
jets have been selected according to their shape as broad jets (black squares ) and 
narrow jets (white squares). The lower part of (a) shows the ratio between the mea- 
sured da/drf^^ for the narrow- and broad-jet samples (triangles). The calculations 
of CDM (solid lines) and MEPS (dashed lines) are included in (a). In (b), the 
calculations of CDM for gluon- (dashed line) and quark-initiated (dot-dashed line) 
jets are included. Other details are as in the caption to Fig. C3 
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Figure 15: (a) Measured differential ep cross-section da/dE^^ for inclusive jet 
photoproduction in the range —1 < rj^^^ < 2.5 in the kinematic region defined by 
< 1 GeV^ and 142 < W^p < 293 GeV. (h) Measured differential ep cross-section 
da/dE^r^^ for inclusive jet DIS in the range —1 < rj^^^ < 2.5 in the kinematic region 
defined by > 125 GeV^. In (a) and (b), the jets have been selected according 
to their shape. In (c) and (d), the photoproduced jets have been selected according 
to the subjet multiplicity and the combination of jet shape and subjet multiplicity, 
respectively. Other details are as in the captions to Figs. (731 and\14\ 
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Figure 16: Measured differential ep cross sections for dijet photoproduction with 
E^^'^ > 17 GeV, E}^'^ > 14 GeV and -1 < rj^^^ < 2.5 in the kinematic region defined 
by Q'^ < 1 GeV^ and 142 < W^p < 293 GeV as a function of (a) | cos6'*| and (h) 
^^^Ko-^d for MJj > 52 GeV, (c) M^J for |cos^*| < 0.8 and (d) The cross 

sections are for events with broad-broad (dots), narrow-narrow (open circles) and 
broad-narrow (stars) dijet configurations selected according to their shape. Other 
details are as in the caption to Fig. El 
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Figure 17: (a) The as{Mz) values determined from the Q CD fit of the measured 
integrated jet shape {ipir = 0.5)) in the different E^^^ regions (squares), (h) The 
combined value of as{Mz) obtained using all the E^^^ regions (square). In both plots, 
the inner error bars represent the statistical uncertainties of the data. The outer 
error bars show the statistical and systematic uncertainties added in quadrature. 
The dotted vertical bars represent the theoretical uncertainties. 
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